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Intrinsic and substrate induced spin-orbit interaction in chirally stacked trilayer 

graphene 

Andor Kormanyo^ and Guido Burkard 

Department of Physics, University of Konstanz, D-78464 Konstanz, Germany 

We present a combined group-theoretical and tight-binding approach to calculate the intrinsic 
spin-orbit coupling (SOC) in ABC stacked trilayer graphene. We find that compared to monolayer 
grapheneSS', a larger set of d orbitals (in particular the (^^2 orbital) needs to be taken into account. 
We also consider the intrinsic SOC in bilayer graphene, because the comparison between our tight- 
binding bilayer results and the density functional computations of Ref. |40| allows us to estimate the 
values of the trilayer SOC parameters as well. We also discuss the situation when a substrate or 
adatoms induce strong SOC in only one of the layers of bilayer or ABC trilayer graphene. Both for 
the case of intrinsic and externally induced SOC we derive effective Hamiltonians which describe 
the low-energy spin-orbit physics. We find that at the K point of the Brillouin zone the effect 
of Bychkov-Rashba type SOC is suppressed in bilayer and ABC trilayer graphene compared to 
monolayer graphene. 

PACS numbers: 73.22.Pr,71.70.Ej,75.70.Tj 
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I. INTRODUCTION 

The low-energy properties of multilayer graphenei de- 
pend crucially on the stacking order of the constituent 
graphene layers^""—. In the case of trilayer graphene, 
there are two stable stacking orders: (i) ABA or Bernard 
stacking and (ii) ABC or chiral stacking. Recent ad- 
vances in sample fabrication methods have resulted in 
high-quality trilayer samples which can be used to probe 
many of the theoretical predictionsi^rJ^. ABC stacked 
trilayer graphene appears to be particularly exciting be- 
cause it is expected to host a wealth of interesting phe- 
nomena, such as chiral quasiparticles with Berry phase 
Stt— , a Lifshitz transition of electronic bands due to trigo- 
nal warpin g^d^ , band-gap opening in an external electric 
ficld^iiii^iSiiiSiiiii^ii^r— , aud broken symmetry phases at 
low electron densities^ir— , to name a few. 

Although there are a number of 
theoretical2rIiiiai2ai22i2^ and experimentali^""ili2^""i22 

studies on the electronic properties of ABA and ABC 
stacked trilayer graphene, the spin-orbit coupling (SOC) 
in these systems has received much less attention. ABA 
trilayer graphene was considered in Ref. [29| within a 
framework of an effective low-energy theory, whereas the 
case ABC stacking was only briefly mentioned in Ref. Isol. 
The understanding of spin-orbit interaction would be 
important to study other interesting and experimen- 
tally relevant phenomena such as spin relaxation^"—, 
weak-localization22i or even spin-Hall effeclj^ in trilayer 
graphene. The recent report of Ref. [isl on the fabrication 
of high mobility double gated ABC trilayer graphene 
may open very promising new avenues for trilayer 
graphene spintronics as well, similarly to the monolayer 
case where highly efflcient spin transport has recently 
been reported^, but with the additional advantage that 
external gates can open a band gap in ABC trilayer 
graphene. 



detailed understanding of the spin-orbit coupling (SOC) 
in chirally stacked trilayer graphene. We start by investi- 
gating the case when the system has inversion symmetry, 
i.e. in the absence of external electric fields, adatoms or a 
substrate. This is the case of intrinsic SOC. The intrin- 
sic SOC opens a band-gap at the band-degeneracy points 
without introducing spin polarization. Previous ah ini- 
tio calculations on monolayer-SIi2S and bilayer— graphene 
provided a strong evidence that the key to the under- 
standing the SOC in flat graphene systems is to take 
into account the (nominally unoccupied) d orbitals in the 
description of electronic bands. Here we take the same 
view and by generalizing the work of Ref. [s^ derive the 
intrinsic SOC Hamiltonian for ABC trilayer graphene. It 
turns out that the most important d orbitals to take into 
account are the dxz, dyz and dz2 orbitals. While the for- 
mer two have been considered in Ref. [s^ in the context of 
monolayer graphene, the latter one is important to un- 
derstand the SOC in AB and ABC graphene. We obtain 
explicit expressions for the SOC constants in terms of 
Slater-Koster— hopping parameters. We also re-derive 
the intrinsic SO Hamiltonian for bilayer graphene^^iiS. 
Through the comparison of our trilayer and bilayer an- 
alytical results with the recent ab initio calculations of 
Ref. 113 we are able to make predictions for the actual 
values of the SOC parameters in ABC trilayer graphene. 
The theory involves electronic bands which are far from 
the Fermi energy Ep but are coupled to the physically im- 
portant low-energy bands close to Ep and hence compli- 
cate the description of the electronic properties. There- 
fore, we derive an effective low-energy Hamiltonian and 
calculate its spectrum. This helps us to understand how 
SOC lifts certain degeneracies of the electronic bands. 

Generally speaking, due to the low atomic number 
of carbon, the intrinsic SOC in single and multilayer 
graphene is weak (according to our prediction, the SOC 
parameters are of the order of 10 /ieV in ABC trilayer, 
the same as in mon» 
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and bilayer---' graphene). 

In this paper we aim to make the first steps towards a Recently however, there have been exciting theoretical 
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proposals to enhance the strength of SOC in monolayer 
graphene and hence, e.g., make the quantum Spin Hall 
stated observable. These proposals suggest deposition 
of indium or thallium atoms^"^ or to bring graphene into 
proximity vifith topological insulators^^. Indeed, very re- 
cently the combined experimental and theoretical work 
ofRef.1551 has provided evidence of a large (10—100 meV) 
spin-orbit gap in monolayer graphene on nickel substrate 
with intercalated gold atoms. Motivated by these studies 
we also discuss what might be a minimal model to de- 
scribe the case where the the SOC is strongly enhanced 
in only one of the layers of bilayer and ABC trilayer 
graphene. 

Our work is organized as follows. In Section |ll] we 
present the tight-binding (TB) model of ABC-stacked 
graphene and introduce certain notations that we will be 
using in subsequent sections. In Section IIIH employing 
group-theoretical considerations and the Slater-Koster— 
(SK) parametrization of transfer integrals, we derive the 
SOC Hamiltonian in atomistic approximation at the K 
point of the Brillouin zone. We repeat this calculation 
for bilayer graphene in Section |IV] so that in Section IVl 
we can make predictions for the actual values of the SOC 
parameters. Next, in Section IVH using k- p theory and 
the Schrieffer- Wolff transformation^^'^^, we derive an ef- 
fective low-energy SOC Hamiltonian which is valid for 
wave-vectors around the K {K') point. Finally, in Sec- 
tion IVHl we consider the case when SOC is enhanced in 
one of the graphene layers with respect to the other(s). 



II. TIGHT-BINDING MODEL 

The basic electronic properties of ABC trilayer are well 
captured by the effective mass model which is derived as- 
suming one Pz type atomic orbital per carbon atom. This 
model has been discussed in details in Refs. |^ and ITTI . 
therefore we give only a very short introduction here (see 
also Fig. [ij . There are six carbon atoms in the unit cell 
of ABC trilayer graphene, usually denoted by Al, Bl, 
A2, B2, A3, B3, where A and B denote the sublattices 
and j = 1,2,3 is the layer index. The parameters ap- 
pearing in the effective model are: 70 for the intra-layer 
Aj — Bj nearest-neighbour hopping, 71 for the interlayer 
hopping between sites Bl — A2 and B2 — A3, 73 (74) de- 
scribes weaker nearest-layer hopping between atoms be- 
longing to different (the same) sublattice, and finally 72 
denotes the direct hopping between sites Al and B3 that 
lie on the same vertical linein the outer layers 1 and 3. 
These hoppings can be obtained by e.g. fitting the nu- 
merically calculated band structure with a TB model^i^. 
The six pz orbitals in the unit cell give rise to six elec- 
tronic bands^iii. As shown in Figd^d), at the K {K') 
point of the Brillouin zone (BZ) two of these bands lie 
close to the Fermi energy E-p — and we will refer to 
them as "low-energy" states. In addition, there are four 
"split-off" states far from Ep at energies E = ±71. 

To obtain the intrinsic SOC Hamiltonian of ABC tri- 




FIG. 1. Lattice and band structure of ABC trilayer graphene. 
a) Lattice structure, where atoms on different layers are in- 
dicated with different symbols, ai and sl2 are the two lattice 
vectors, b) Schematic side view of the unit cell with the most 
important hopping amplitudes, c) Schematic of the Brillouin 
zone with reciprocal lattice vectors bi, b2, the high symme- 
try points r, K and K' and the momentum p = (Px,Py) 
measured from K. d) Schematic of the band structure at the 
K point of the Brillouin zone. The velocity «o is given by 
vo — -\/3/2a7o/ft, where a is the lattice constant: a = |ai|. 

layer graphene we generalize the main idea of Ref. [s^ 
where the SOC of monolayer graphene was discussed. 
Using group-theoretical considerations and density func- 
tional theory (DFT) calculations it was shown in Ref. [s^ 
that in the case of monolayer graphene the main contri- 
bution to the intrinsic SOC comes from the admixture of 
Pz orbitals with some of the (nominally unoccupied) d or- 
bitals, namely, with the dxz and dyz orbitals. The other 
d orbitals, dz'^, dxy and dj.2_j,2 play no role due to the 
fact that they are symmetric with respect to the mirror 
reflection cr^ : {x,y,z) (x,y,—z) to the plane of the 
graphene layer, whereas Pz is antisymmetric. The impor- 
tance of the d orbitals for the understanding of SOC in 
monolayer graphene was also pointed out in Ref. [H. 

The symmetry group of ABC stacked trilayer graphene 
is R3m {D^d) which does not contain the mirror reflection 
Uh- Therefore, from a symmetry point of view, in an 
approach similar to Ref. [H, all d orbitals need to be 
taken into account. To derive a TB model one should 
therefore use as basis set the Bloch wavefunctions 

*;^(r,k) = ^5]e'''-(^"+*")$^^(r- (R„+t„)), (1) 

where the wavevector k is measured from the F point 
of the BZ (see Fig. [T|), aj = {Al,A2,A3,Bl,B2,B3} 
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is a composite index for the sublattice a = {A, B} and 
layer j = {1,2,3} indices and denotes the atomic or- 
bitals of type /3 = {p^, ^^2, doj^, dyz,dxy, d^2_y2} in layer 
j. The summation runs over all Bravais lattice vectors 
R„, whereas the vectors t„ give the position of atom a in 
the two-dimensional unit cell. We use a coordinate sys- 
tem where the primive lattice vectors are ai — |(\/3, 1) 
and a.2 = ^{—^/3,1), the positions of the atoms in the 

unit cell are t^i = tss = (0,0), tA2 = t^i = § 

and t^3 = ±32 = § ("73' l) where a — 2. 46 A is the lat- 
tice constant. The K and K' points of the Brillouin 
zone, which are important for the low energy physics 
discussed in this paper, can be found at K = (0, |^), 

The symmetry group of the lattice contains three- 
fold rotations by ±27r/3, about an axis perpendicu- 
lar to the graphene layers. Since the atomic orbitals 
dxz,dyz,dxy,dx2^y2 thcmselvcs do not possess this sym- 
metry, instead of VE'^-' (r, k) given in Eq. ([1]) we will use 
Bloch states which depend on -^{d^z ± idyz) ^ ^Y^^, 
-^{dx2-y2 ± idxy) ^ Y^'^ (rotating orbitals), where 
Yj™ are spherical harmonics. Taking into account that 
Pz ~ and dz2 ~ Y2 , the Bloch states we use as basis 
will be denoted by ^'"toj (r, k), where Z = 1, 2 and mi = 0, 
whereas 1712 can take all allowed values m2 — —2 ... 2. 
Often, we will need a linear combination of two of these 
basis functions where both of the basis functions have 
the same quantum number / but one of them is cen- 
tered on an A type atom and the other one is on a 
B type atom, e.g. ^[*f_i(r,k) - «'fi](r,k)]. As a 
shorthand notation, we will denote the symmetric com- 
bination of two such basis functions by ^'f ;^ (r, k) = 
--^[\['^j^(r, k) -I- 4'^^,(r,k)] and the antisymmetric one 
with ^i:i,^,(r,k) = ^[vl';^/„(r,k) - *f,{,(r,k)]. The 

first upper index j in 5';';^^ always denotes the layer 
index of the atomic orbital centered on the A type atom, 
the second upper index j' is the layer index for the or- 
bital centered on the B type atom, the first lower index 
I is the common angular momentum quantum number, 
finally, the second and the third lower indices m, to' give 
the magnetic quantum number in the same manner as 
the upper indices give the layer index. To lighten the 
notation, we will usually suppress the dependence of the 
Bloch functions on (r, k) and use the bra-ket notation, 

e-g- l^lio): l^>- 

The derivation of the spin-orbit Hamiltonian pro- 
ceeds in the same spirit as in Ref. H^: (i) first, we 
neglect the spin degree of freedom. Using the Slater- 
Koster parametrization to describe the hopping integrals 

(*lm,m'l^l*;'','m'''m"') ("^ the single particle Hamil- 
tonian of the system) at a high symmetry point (the K 
point) of the Brillouin zone and group theoretical consid- 
erations we obtain certain effective Bloch wavefunctions 



which comprise Pz and d atomic orbitals centered on dif- 
ferent atoms; (ii) using these effective wavefunctions we 
calculate the matrix elements of the spin-orbit Hamilto- 
nian in atomic approximation, and (iii) employing the 
k- p theory we obtain the bands around the K point and 
then we derive an effective low energy Hamiltonian. 



III. INTRINSIC SOC 

If, in addition to the Pz orbitals, we include also the d 
orbitals into our basis, there will be six basis functions 
I^Tm) centered on each of the six carbon atoms in the 
unit cell and hence the TB Hamiltonian HabCj which 
is straightforward to calculate in the SK parametriza- 
tion, is a 36 X 36 matrix. By e.g. numerically diago- 
nalizing this matrix one would find that the Pz orbitals 
hybridize with some of the d orbitals and one could see 
how the low energy and the split-off states, obtained in 
the first instance by neglecting the d orbitals, are modi- 
fied. According to band theory each state at the K point 
should belong to one of the irreducible representations 
of the small group of the K point^, which is 32 (D3) 
in this case. This group has two one-dimensional irre- 
ducible representation, denoted by Tai and F^a respec- 
tively, and a two-dimensional one denoted by F^; (see 
Appendix The matrix elements of Habc between 
basis states corresponding to different irreducible repre- 
sentations of are zero—-. In other words, Habc can 
be block-diagonalized by choosing suitable linear combi- 
nations of the basis functions such that the new 
basis functions transform as the irreducible representa- 
tions of the group D3 because the hybridization between 
Pz and d orbitals will preserve the symmetry properties. 
A group-theoretical analysis of the problem shows that 
in a suitable basis Habc is block-diagonal having (i) two 
6x6 blocks which we denote by Hr^^.^ and -ffr^a ; they 
correspond to basis states with F^^ and Ta2 symmetry, 
and (ii) there is one 24 x 24 block Hr^ corresponding to 
states with F e symmetry. (The basis vectors with F^^ , 
Ta2 a-iid Te symmetries are listed in Appendix [AI Table 
IVIip . As a concrete example we will consider i?rAi 
discuss how one can extract an effective orbital in which 
Pz atomic orbitals with large weight and d orbitals with 
small weight are admixed. The calculation for -ffr^a ^^'^ 
Hye cases is analogous and will be presented only briefly. 

The basis states transforming as the irreducible repre- 

sentation F^^ are l*};;^ ^), 1^-^;} \^^2lfi). l*2;-2,2), 

l*2;-i,i) and |^'2;2 -2)- The TB Hamiltonian i/r^, can 
be further divided into 3x3 blocks: 

/ H^} W^^ \ 
Explicitly, the upper left block H^j , corresponding to 
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the basis states |^'};o,o)' l*2!i -i) and I'^l'^^o), reads 



Pz-hke state I^'x'q q) is transformed into 



H. 



( £0 — 72 -^VpdTT 



pd 



A1,B3 



_3_T/ 

^ ypdiT 

1 rrAl,B3 

\ ^pda 



£d 








(3) 



A1,B3 
dda 



where the upper indices Al^ _B3 on the SK parameters in- 
dicate the atomic sites between which the hopping takes 
place. The parameter Vpd-w describes hopping between A 
and B type atoms within the same graphene layer and 
we assume that its value is the same in all three lay- 
ers. The matrix elements in W^^ of Hamiltonian ([2]) are 
either zero^ or describe skew hoppings between the pz 
and d orbitals located on different atoms. We assume 
that these skew hoppings are much smaller than both 
the vertical hopping V^^^^ and Vpd-n ■ This is not a cru- 
cial assumption and the neglected skew-hoppings can be 
taken into account in a straightforward manner. How- 
ever, it simplifies the lengthy algebra that follows and 
we believe it yields qualitatively correct results (see Sec- 
tionlIV[). With W^'^ « we see that H^j (corresponding 

to the basis functions |*2i-2,2); l^2!-i,i) and {^'2'^, -2)) is 
decoupled from Hpd and that this latter matrix describes 
the hybridization between the pz orbital based basis vec- 



tor l^'i'oo) and the basis vectors 1^*2'^ 
volving dxz, dyz and dz2 orbitals. By diagonalizing H;^J^ 
one could find out how the energy = Ep — 72 of one 
of the low-energy states is modified by the d orbitals. 
The secular equation det{H^J^ — el) = leads to a cu- 
bic equation in e and the solutions can only be expressed 
using the Cardano formula. Instead, we next perform 
a Schrieffer- Wolff transformation (Lowdin partitioning) 



3,1 



a), 1*2:0,0) 



H e to approximately block-diagonalize 



^pd - 

^pd i^^*^ a 1 X 1 and a 2 x 2 block by eliminating the ma- 
trix elements between \'^\'q q) on one hand and |^'2'i 
1*2 0) '^'^ other hand. (A detailed discussion of this 
method can be found in e.g. Refs . l46l and l47l. ) The matrix 
S is anti-Hermitian: — —S and only its non-diagonal 
blocks Spd and —S^^ are non-zero. In first order"*^ of the 
coupling matrix elements Vpd-R and V^^^^^ one finds that 



Vr 



pdn 



A1,B3 
pda 



Se 



pd 



72 Sspd - 72 



(4) 



where Vpd^r = 75 ^pd^' ^^pd 



Ep - Ed and 72 = 72 



A1,B3 
dda 



The 1 X 1 block of reads EA^ 



12+V^d^/ {Sspd - y2) + ivf^];^^y/iSEpd - 72), this means 
that the energy ea^ of the low energy state is shifted 

by VX/iSspd - 72) + (V^^'^'r/iSEpd - 72). While the 
Schrieffer- Wolff transformation is usually used to obtain 
effective Hamiltonians, one can also obtain the new basis 
in which H;^J^ is block-diagonal. Making the approxima- 
tion e""^ « 1 — 5 (see Ref. [iol ) we find that the purely 



I* 



1,3 X 
1,0,0/ 



Vpd-n 



Se 



pd - 72 



I* 



A1,S3 
pdcr 



2,1,-1/ 



5e- 



pd - 72 



I* 



1,3 
2,0,0 



(5) 

i.e. it is admixed with two other basis vectors contain- 
ing dxzidyz and dz^ orbitals. Since l^'p^ ) corresponds 

to a 1 X 1 and hence diagonal block of Hpd, it is an 
approximate eigenvector of H^J[ with energy ea^ ■ The 
upper index pz in l^'p^ ) is meant to indicate that in 
this state Pz orbitals have the largest weight. There are 
two other states with T Ai symmetry which could be ob- 
tained by diagonalizing the the remaining 2x2 block of 
H^^ . In these states |^'2:i,-i) and |^2io,o) would have 
large weight. They are however far remote in energy 
from l^'p^ ) and therefore play no role in our further 
considerations. The situation will be similar in the case 
of the two other irreducible representations, T A2 and Fg, 
therefore we will suppress the upper index pz henceforth 
in the notation of the physically important approximate 
eigenstates. 

We now briefiy discuss the symmetry classes T A2 and 
r^;. The calculation for the other 6x6 block of Habc 
with Ta2 symmetry is completely analogous to the Tai 
case, the resulting approximate eigenvector, l^'p^ ) is 
shown in the left column of Table IIIII Its energy, apart 
from the shift due to the d orbitals, which will be ne- 
glected, is EA2 — £p + l2- 

The matrix block corresponding to states with T e sym- 
metry can be written as 



Hi 



W^d 



TT I ' PP 



dd 



(6) 



Here the 4x4 block H^^ contains the matrix elements be- 
tween the basis vectors \-^l'lo)^ l*?;o,o)> l*?;o,o)> l*?;o,o): 

the 4 x 20 block is the coupling matrix between the 
above shown pz orbital based basis vectors and the d 
orbitals based basis vectors. (The full set of basis vec- 
tors for each of the three symmetry classes is listed in 
Appendix [X] ) A direct calculation shows that Hpp is a 
diagonal matrix with Ep -f 71 , ep -f 71 , ep — 71 , ep — 71 
entries. Earlier we have referred to these states as "split- 
off" states. Since H^^ is diagonal, its approximate eigen- 



vectors hl'T 



and |*r 



which 



are listed in Table IIIII can be obtained in exactly the 
same way as I^Ta )■ 

We will refer to the basis formed from the 
physically important approximate eigenvectors 



E2.1 



)J*r._)} as 



"symmetry basis" henceforth. (The symmetry basis for 
the K' point can be obtained by complex-conjugation.) 
Looking at these basis vectors we see that in contrast to 
monolayer graphenc'^^, where only dxz and dyz orbitals 
hybridize with the Pz orbital, here also the dz^ orbitals 
are admixed. As it will be shown below, the admixture 
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of dz2 orbitals is crucial to obtain the non-diagonal 
elements of the SOC Hamiltonian. 


SOC 


*c1f 










*c1l 


We can now proceed to calculate the SOC Hamilto- 
nian. This can be done in the atomic approximation, 




Al/2 Sz 





A2 s+ 


Ais; 


A4 SI 


A3 SI 


whereby the spin-orbit interaction is described by the 
Hamiltonian 







— Ai/2 Sz 


-A3S; 


— A4 s+ 


-AiSI 


-A2 SI 


Tjatomic ^ ^ dV (r) 

~ Amlc? r dr ^ ' 


>T|S1 


\ 

Ai SI 


—A3 6_ 
-A4 SI 


— A3/3 ^2 




u 

A|/3 Sz 


— ^3/3 J + 




u 

A3/3 s+ 


Here V{r) is the spherically symmetric atomic potential, 
L is the angular momentum operator and S = {Sx, Sy) is 




A4 


-Ais; 


^A3/3 SI 





^'^3/3 'S'z 





a vector of spin Pauli matrices Sx^ Sy (with eigenvalues 




A3S; 







A3/3 SI 





A3/3 Sz 



±1). Introducing the spinful symmetry basis functions by 
l^*^) |5'^(8)s), where s = {'[, \.} denotes the spin degree 
of freedom and noting that L- S = LzSz + L^S- + L^S+, 
where L± = Lx±iLy and S± = ^{Sx^iSy), it is straight- 
forward to calculate the matrix elements (iJ^^'^)^^^ — 
in the symmetry basis introduced ear- 
lier. Using the notation 5'J = ■^{Sx ± ir- Sy), where 
T = +1 (—1) corresponds to the K [K') point of the BZ, 
the result is shown in Table HI 



SOC 



^1/30+ ^1/3''- ^-^1/3''+ -^1/3 



Al/2Sz 

\ El QT 

'^1/3'^- 

\ El QT \ El QT ri 

'^1/3'^+ -^2/3''+ " 



^2/3''+ '^2/3'^- ^2/3''+ ^2/3'^- 



\El QT 

'^2/3''- 



'2/3 






2/3*- 



'2/3^ 



\E2 QT \E2 QT \z Q \ or 

"'^1/3'-'- ^2/3"^- ^3/3 ^3/3J + 

\ ^2 QT \ E2 QT \ QT \Z Q 

'^1/3'^+ -^2/3''+ ^3/3'5- ^3/3'-'z 



A|_/3Sz A3/3S+ 

A3/3SI A|y3Sz 







TABLE I. Intrinsic spin-orbit Hamiltonian r Hf^'~' in the 
symmetry basis. Here r = 1) corresponds to the K 

(A'') point. 



TABLE II. tH^J^'^ in the basis of effective pz orbitals. 



The physical meaning of these SOC parameters is 
probably more transparent if one rotates Hf^'^ into the 
basis of on-site effective Pz orbitals (see Appendix IB ip 
. These basis vectors result from the admixture of an 
on-site pz orbital |5'"o) with large weight and l^'j ±1), 
jvl/^i)) (a 7^ a') with small weight. H^^^^ in the on-site 
effective pz basis is shown in Table HIl The SOC param- 
eters Ai,A2,A3 and A4 are linear combinations of A^y^g, 

A03, A^^ and A^^^^, the explicit relations are given in 
Section |Vl In Section El we will also make predictions 
which might be useful to guide the fitting procedure if 
results of DFT calculations are fitted with a TB model, 
as e.g. in Ref. EO. 

Since the lattice of Bernard stacked bilayer graphene 
has the same symmetry group as ABC trilayer, the con- 
siderations made in this sections can be easily applied to 
bilayer graphene as well. A brief summary of the bilayer 
calculations is given in Section IIVI The importance of 
the bilayer results is that they can be compared with the 
numerical calculations of Ref. H^- Based on this compar- 
ison we will be able to estimate the values of five of the 
seven SOC parameters of ABC trilayer. 



The SOC Hamiltonian shown in Table U is the main 
result of this Section. Explicit expressions in terms of 
SK parameters for the coupling constants appearing in 
Table [T] can be found in Table IIIII In contrast to previ- 
ous works where SOC in ABC trilayer was discusse d"^"'^^ , 
we find that the number of SOC parameters is seven^. 
As we will show, , X^J^ , and X^J^ are related to 
inter-layer SOC and calculations which are based on the 
symmetry properties of low-energy effective Hamiltoni- 
ans may not capture them. The A3/3 parameter ensures 
that the otherwise four-fold degeneracy of the split-off 
states at the K point is lifted, as it is dictated by general 
group-theoretical considerations^^ii^. These five param- 
eters are proportional to the product Vpd-nVpda and they 
could not be obtained considering only the dxzi dyz or- 
bitals and in-plane SOC. The remaining two SOC param- 
eters, Ai/2 and Ag^g are proportional to V^^^ and describe 
in-plane SOC. 



IV. INTRINSIC SOC IN BILAYER GRAPHENE 

In this Section we give a brief summary of our TB 
calculations for the intrinsic SOC in bilayer graphene 
and compare the results to the DFT computations of 
Ref. The low-energy states of bilayer are also found 
at the K and K' points of the BZ, hence the cal- 
culation follows the same steps as in Section IIIII (i) 
first we obtain the basis states of the symmetry basis 
{|*r^,),|^'r^,),|*rB,,,),|*rB^ ^)}, ^nd (ii) we calculate 

the matrix elements of 7?^^™"^. Note that in the case 
of bilayer graphene there is only one pair of bands which 
transforms as the two-dimensional representation Te- 
For easier comparison we adopt the notation of Ref. 
for the SOC parameters. As it is shown in Table ITVl the 
SO Hamiltonian in the basis of the effective Pz orbitals 
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Al/2 

^1 
/3 



^1/3 



2/3 
^2/3 



■^3/3 



9 f 

2 '''' (<Sepd+72){<Sepd-72) 



A3/3 = 3\/3^c 



(*epd+7l) 


\ ^epd-72 




iepd-72 


^pd. 


/ Vpd„ 


+ 




(*Epd 11) 


\ '5epd-72 


*£pd-72 




/ %da 


+ 


^pdcr 


(*epd+7l) 


\ '5epd+72 


'5epd+72 


Vpd„ 


/ V-pd. 




^pd<T 


(*epd-7l) 


V '5epd + 72 




iepd+72 




^pdCT 







('5Epd+7l)(iepd-7l) 



9 c {Vpd^r 
'2 



(*epd+7l){*epd-7l) 



TABLE III. Symmetry basis functions (left column) and intrinsic spin-orbit matrix elements (left column) in terms of Slater- 



Koster parameters. Here is the angular momentum resolved atomic SOC strength, Vdda = V^, 



Vpda = V 



T/ _ ■\/A2,B\ _ ■,rA3,B2 
yvda y^rl^ pda ' 



Sspd = Ep - Ed, 72 = 72 + Vddcy and we assumed that Vpd-^ = V > 



= 1, 2, 3. Although the basis 

functions shown in the right are not normalized, the SOC parameters are correct in the lowest order of the products of the 
small parameters Ypda/iSspd ± 72), Ypda/iSspd ± 72), Ypd^v/ {&£pd ± 71) 



can be written in a form which, apart from a unitary 
transformation, agrees with the result given in Table IV 
of Ref . [13 (see Appendix |X] for details) . In terms of the 



SOC 



Ao 5+ 
SI 




-All Sz 


\ bi QT 







\ hi QT 
A4 ^ + 








-Ai' SI 


An Sz 


-\oSL 


— Ao S+ 


— A12 Sz 



TABLE IV. Intrinsic spin-orbit Hamiltonian r of bilayer 
graphene in the basis of effective Pz orbitals. 



A4 = — 12/ieV. These numbers are in qualitative agree- 
ment with the predictions that we made below Eq. 
for the SOC parameters. If, in addition, one assumes 
that VpdTrVpda < then according to ^ three of the pa- 
rameters (All, A12, Ao) should have the same sign, which 
would again agree with the results of Ref. [iOi 



SK hoppings, the SOC parameters read: 

VpdTrVpda 71 



Ao = 3V3^d 



An 



9 ^ ^pdiT 



pd 



\bi _ o fn c ^pdTrVpdn _ 
A4 — •3VO?dTT2 

"^pd 



~ 2 ' 

71 



A 



12 



9 



(8) 



where 71 = 71 -I- Vdda ■ Looking at the expressions given 
in ([5]), one can make the following observations: (i) since 
Sspd = Ep ~ Ed < and 71 > 0, one can expect that 
the sign of Aq and A4' will be different; (ii) |Ao| < |A4*| 
because < 1; (iii) An and A12 have the same sign and 

|Cpd I 

they are approximately of the same magnitude because 
7iVfepd < 1- 

By fitting the band structure obtained from DFT cal- 
culations with their tight binding model, the authors of 
Ref. HO found the following values for the bilayer SOC 
parameters: 2Aii = 24/ieV, 2Ai2 = 20/xeV, Aq = 5fieV, 



V. SOC PARAMETERS FOR TRILAYER 
GRAPHENE IN TERMS OF SK HOPPINGS 



We are now ready to make predictions for five of the 
seven ABC trilayer SOC parameters. To this end, we 
first express the SOC parameters in the effective Pz or- 
bital basis in terms of the SOC parameters obtained in 
the symmetry basis. Moreover, using the formulae given 
in Table IIIII for Xfj'^ , A^Jg , ^ys ^^'^ ^2/3 ' '^^^ '^^^ ^^^'^ 
express Ai, A2, A3 and A4 in terms of the SK hoppings. 
Taking into account the results of the previous section, 
this then facilitates to make predictions for the numerical 
values of the ABC trilayer SOC parameters. 

First, the SOC parameters in terms of the SK hop- 
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pmgs: 



Ai 



= (■^1/3 + -^2/3 ~ ^f/3 ~ ^f/3)/2 
72 



-Vpda + Vpdij-?-— 

oepdj 



1] 

/3 + ^^1/3 ^2/3 '^2/3''/^ 



^2 - (Af/3 

~ n \ VpdaY^ + VpdaY^ ] . 
V OSpd 0£pdj 

^ + -^f/3 + ^2/3 + ^2/3)/^ 

71 72 
7 H V'pdo- 



V OSpd 0£pd J 

A4 = (Af/3 + Af/3 + Af^^3 + Af/3)/2 



. , ,7 71 ^ 

' '7 ^pdcrT 1- Vpda-. , 



(9a) 



(9b) 



3c) 



(9d) 



where 77 = 3-^3 



Similarly to the bilayer case, 



looking at Table Hill and Eqs. ^ one can make the follow- 
ing observations: (i) one would expect that A1/2 ~ A3/3, 

(ii) IA3/3I « |Ai| and assuming that Vpda/Vpdo oc 71/72 
one finds that |Ai| > IA2I > IA3I > IA4I, (iii) A4 has op- 
posite sign than A3 and similarly for A3/3 and Ai because 
the second term in the expression for Ai in Eq. (j9ap can 
be neglected with respect to the first one. 

Comparing the expressions in terms of the SK hoppings 
given in Eqs. ([9]) with the corresponding ones for the 
bilayer case in ([5]), the following estimates can be made: 



2A 



1/2 



2AI 



3/3 ^ 20/ieV, IA3/3I w I All w lO^eV and 
IA2I ~ 5/ieV. Since A3 and A4 are proportional to Vpda 
(assuming Vpda/ypda 71/72) which is unknown, we 
cannot give a numerical estimate for their value. One 
would expect that they are much smaller than Ai and A2 
because Vpda corresponds to a remote, and presumably 
weak p — d hopping between the Al and B2> sites. 



VI. EFFECTIVE SOC HAMILTONIAN 

The calculations in the previous sections are valid, 
strictly speaking, only at the K point of the Brillouin 
zone. To obtain the Hamiltonian in the vicinity of the 
K point, where the states close to the Fermi energy can 
be found, one can perform a k- p expansion of the bands. 
We neglect the weak k dependence of the SOG^, hence 
the total Hamiltonian of the system can be written as 
Habc - H^.l"^ + tH^J^^. Here H^,^^ is the k- p 
Hamiltonian obtained without taking into account the 
SOC, whereas Hf'^^ is the SOC Hamiltonian calculated 
at the if point. H^,^^ has been published before, see 
e.g. Refs. 13 and [ill 

To study the low energy physics however, in which we 
are primarily interested, the use of H^^^ is not con- 
venient, since it includes four bands that are split-off 
from the Fermi energy of the (undoped) ABC trilayer 



by the large energy scale ~ ±71— ill. Therefore we de- 
rive an effective two component (or, including the spin, 
four component) Hamiltonian H'^^q which describes the 
hopping between atomic sites Al and BZ. To this end 
we again employ the Schrieffer- Wolff transformation and 
keep all terms which are third order or less in the mo- 
mentum TT, tt''' and first order in the SOC constants. Here 
TT = —{ipx + Tpy), where r = 1(— 1) for valley K [K') 
and the momenta Px^Py are measured from the K {K') 
point of the BZ, see Fig. [TJc). Keeping terms up to third 
order in tt, tt^ is essential to reproduce the important 
features of the low-energy band structure^iiii, such as the 
band degeneracy and the trigonal warping. The neces- 
sary formulae for the matrix elements of the effective 
Hamiltonian can be found in Ref. [4^ We treat 71 as 
a large energy scale with respect to v^-n, I73I, 74, 72 and 
the typical energies E we are interested in and keep only 
the leading order for the terms involving 72, and 114. 
(The velocities Vi are given by Vi — {^/3/2)a'-fi/h, where 
a — 0.246 nm is the lattice constant of graphene.) For the 
folding down of the full Hamiltonian we use the form of 
Habc hi the symmetry basis because in this case all the 
large matrix elements are on the diagonal and therefore 
the quasi-degenerate perturbation approach is expected 
to work well. Once we obtain the effective Hamiltonian 
H'abc ill symmetry basis we rotate it into the basis 
of effective pz orbitals centered on atomic sites Al and 
53 because assumes a simpler form in this basis. 
Explicitly, one can write H'abc ~ ^tf + , where the 
electronic part is given by 



Hchi 



7? I 7r3 



VOV3 r t , tl ( 1 

71 M 1 



Hj2 = 72 1 



Hv 




1 

1 



71 \ nn^ 



(10a) 
(10b) 
, (10c) 
(lOd) 



We note that when applying the Schrieffer- Wolff trans- 
formation we did not assume that tt and tt^ commute, 
therefore the Hamiltonians pH)) and pT|) are valid in the 
presence of finite external magnetic field as well. In zero 
magnetic field, using the notation p^ = tttt^ = tt^tt, the 
Hamiltonian (1101) simplifies to the corresponding results 
in Refs. |8i and 1 111. 
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(11a) 
(lib) 



The effective SOC Hamiltonian is 

rreff „ frmn , rV(l) , fT(2) , A-(3) 
^^so so ^ ^-^ so ^ so ^ so ; 

7i 

rr(2) _ '"0^3 ^ rcr / tN2 I cr 21^ 
-"so -^^^'^3/3['3+(7r') +5_7rjcr2, 



-^i? = -'^ :^(-^3/3 ~ Al/2)'5'z 



Here the Pauli matrix acts in the space of {Al, i?3} 
sites and Ai = Ai + A3(i'3/t;o) + A2(t'4/t'o) « Ai. At 
low energies v^]? ^VQV-j,'p^ ^ 7^ and the corresponding 
terms in pT|) can be neglected. The first term, 77™" is 
the well known SO Hamiltonian of monolayer graphene -^ 
and describes the leading contribution to SOC. The next 
term, H')to is the most important momentum dependent 
contribution close to the K point. Keeping only i^i™"' 
and i/io'' the effective SOC Hamiltonian can be written 




in a more compact form as 



Ai/25. + Ai^(Sxp)^ 
71 



(12) 



where p = {px,Py)- We note that Eq. (fT2|) also describes 
the effective SOC Hamiltonian of bilayer graphene with 
Ai/2 replaced by —An and Ai = — A|* + Ao(w4/wo) 
(for All, Ao and A4* see Section |TV|. 



\ bi 
A4 



In zero external magnetic field H 



off 

ABC 



IS eas- 



ily diagonalizable. Keeping only the leading terms 
in Hgoj we obtain the eigenvalues E± = 
_2B>2l±p'^ ± r(p) (each doubly degenerate) where r(p) = 



71 
•^1/2 



72 — ^^"'^ P + ■:^P e where p = |7r|) whereas 0p is 
the phase of tt. The main effect of SOC on the spectrum 
is, similarly to monolayer^ and bilaye r^°i"^° graphene, to 

open a band gap E^g = 2^X^^^ + |c(pd)p at the band 

degeneracy points pd, while preserving the spin degener- 
acy of the bands. Comparing E^g to the SOC band gap 
in monolayer and bilayer graphene at the K point, we 
expect that in ABC trilayer it should be somewhat big- 
ger due to the |c(pc()p term, i.e. because the band gap 
can be found away from the K point at finite p^. 



VII. SUBSTRATE INDUCED SOC 

The calculations of Section lllll - |V] suggest that the 
intrinsic SOC in ABC trilayer is relatively small, the or- 
der of magnitude of the SOC parameters is lO/ieV, the 
same as in monolayer— "—1^ or bilayer— graphene. One 
can, in principle, enhance the SOC in a number of ways, 
e.g. by impurities^, by applying strong external elec- 
tric fieldffl, making the graphene sheet curved^, using 



adatoms with large atomic number-- or bringing the tri- 
layer into proximity with a suitable substrat o^"^'^^ . Since 
proximity to a substrate or adatoms are likely to lead to 
a larger SOC effect than what one can induce by an ex- 
ternal electric field, here we consider an effective model 
whereby strong SOC is induced in one of the outer lay- 
ers of trilayer graphene whereas SOC is not altered in 
the other two layers. For concreteness, we assume that 
it is the first graphene layer where strong SOC is in- 
duced and for simplicity we will refer both to the sce- 
nario involving a substrate and that involving adatoms 
as "substrate induced SOC". From symmetry point of 
view, the model we consider here is not exact: a group 
theoretical analysis^Sf^i^ of the matrix elements of the 
SOC Hamiltonian shows that by breaking the inversion 
symmetry there can be in principle 21 different SOC pa- 
rameters. However, we assume that all intra and inter- 
layer SOC parameters will remain small with respect to 
the SOC parameters in the layer that is in the immedi- 
ate proximity of the substrate. The relevant part of the 
SOC Hamiltonian that we consider as a minimal model 
is shown in Table |Vl We assume that at the K point of 
the BZ only four SOC parameters have significant values 
and neglect all other intra or inter-layer coupling spin- 
orbit parameters. The SOC parameters that we keep 



SOC 



\l/2Sz \brS\ 
-A1/2S, 



|c(p)P + |d(p)|2, c(p) = -Ai^pe'^P, d(p) = 



TABLE V. The non-zero part of the model SOC Hamiltonian 



are A1/2 and Ag^g which may describe enhanced diagonal 
SOC on AI and Bl type atoms, respectively, whereas 
Xbr is the Bychkov-Rashba^ type SOC acting only in 
the first graphene layer. Since we are going to derive an 
effective low-energy Hamiltonian, we also keep non-zero 
the intrinsic SOC parameter A1/2 on atom Bi, which, 
however, might be much smaller than \if2 and Ag^g. 

Similarly to Section IVTl for the full Hamiltonian of the 
system in the vicinity of the K point we take Habc = 
Htp^ + Htf. In the Hamiltonian H^.'^^ , as an addi- 
tional substrate effect, we also include a possible shift A 
of the on-site energies of atoms AI and Bl with respect 
to atoms in the other two layers. (Otherwise -ffj^p*^ is the 
same as in Section IVll ) To study the low energy physics 
we again use the Schrieffer- Wolff transformation to elim- 
inate the high-energy bands. The hopping 71 is treated 
as a large energy scale with respect to vqit, I73I, 74, 72, 
A and we keep only the leading order for the terms in- 
volving A, 72, ^3 and W4. The electronic part H^f of 
the effective Hamiltonian H'^^fj contains one new term 
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in addition to the terms in Eq. (ITOl) : 



Ha = A 




27? 



2tt^tt + tx'k'^ 

— TTTI 



(13) 

Up to linear order in the momentum and the SOC 
parameters the effective spin-orbit Hamihonian has the 
fohowing terms: 



fjABC 

Aabcaq) 

so 



fjABCiO) , fjABCil) , fTABC,(2) 

' A 



-tA 



BR 




Sl-K 




ffABC,(2) ^\ "^/QTi.QT^ 

71 



It is interesting to compare this resuhs to what one 
obtains for bilayer graphene using the same model^: 



fjAB,{l) , frAB,(2) 




fjABXl) 



71 




(15a) 
(15b) 
(15c) 



In the appropriate limit Hamiltonian (jl5p agrees with 
the result of Ref. HO- Similar models have also been con- 
sidered in two recent works: in Ref. '59' strong Rashba 
SOC was assumed in both layers of bilayer graphene and 
a linear-in-momentum low-energy SOC Hamiltonian was 
derived, whereas in Ref. [s^ the authors kept only the first 

= 0. 



term, H^o^'^^^ and set A: 



3/3 

One can see that in contrast to monolayer graphene, in 
bilayer and ABC trilayer graphene the Rashba SOC af- 
fects the spin-dynamics in the low-energy bands through 
terms which are momentum dependent^S. This means 
that at the K {K') point the effect of Rashba-type SOC 
is suppressed with respect to monolayer graphene. Fur- 
thermore, noting that S> W3, a comparison of the pre- 

factors of Hfo^'^''"^^ and H^o^''^^^ may suggest that the 
influence of linear-in-momentum terms on spin-dynamics 
might be more important in bilayer than in trilayer. This 
is, strictly speaking, only true in the model that we used, 
i.e when other off-diagonal SOC parameters can be ne- 
glected with respect to As^. In general, there would be a 
linear-in-momentum SOC Hamiltonian with a pre-factor 
proportional to 1^0/71 for trilayer case as well. The 
model introduced above could be relevant e.g. in an ex- 
periment similar to Ref. [5^ if bilayer or trilayer is used 



instead of monolayer graphene. Varykhalov et al^ re- 
ported a large Rashba SOC in monolayer graphene with 
Xbr between 10—100 meV, whereas, as we have seen, the 
intrinsic SOC parameters are typically of a few 10/ieV. 

A detailed study of the properties of Hamiltonians 
and ([T5|) is left for a future study. We expect a rich 
physics emerging from the interplay of diagonal and off- 
diagonal spin-orbit terms and the interlayer asymmetry 
A. 



VIII. CONCLUSIONS 

In conclusion, we studied the intrinsic and substrate 
induced spin-orbit interaction in bilayer and ABC tri- 
layer graphene. Assuming that in flat graphene systems 
the most important contribution to the SOC comes from 
the admixture of Pz and d orbitals and using a combi- 
nation of group-theoretical and tight-binding approaches 
we derived the intrinsic SOC Hamiltonian of ABC tri- 
layer graphene. In contrast to the similar calculations 
for monolayer graphene^, we found that in bilayer and 
ABC trilayer in addition to dxz and dyz orbitals also d^a 
orbitals have to be taken into account. For both bilayer 
and trilayer graphene we obtained explicit expressions for 
the SOC parameters in terms of SK hopping parameters. 
By comparing these expressions with the DFT calcula- 
tions of Ref. 140, we were able to estimate the values of the 
intrinsic SOC constants for ABC trilayer graphene. Since 
the intrinsic SOC is quite small, we considered a situation 
when adatoms or a substrate can induce a strong SOC 
(intrinsic diagonal or Rashba type off-diagonal) in only 
one of the layers of bilayer and ABC trilayer graphene. 
To describe the low-energy physics we derived effective 
Hamiltonians for both systems. We found that the effect 
of Rashba type SOC is suppressed close to the K {K') 
point with respect to monolayer graphene. 

The approach that we used here to derive the SOC 
Hamiltonians can be employed in the case of other re- 
lated problems as well. For instance ABA stacked trilayer 
graphene or graphite can be treated on the same foot- 
ing when one takes into account, that they have differ- 
ent symmetries than bilayer and ABC trilayer graphene. 
Considering the substrate induced SOC, which can be 
strong enough to make the SOC related phenomena 
experimentally observable, one interesting question is 
whether the different symmetries and band structure of 
ABC and ABA trilayer would manifest themselves in e.g. 
significantly different spin-transport properties. 
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Appendix A: Small group of the K point and basis 
functions 

The relevant symmetry group for the calculations at 
the K point for both trilayer and bilayer graphene is the 
group 32 (-D3). The generators of this group are two 
three fold rotations (C3) around an axis perpendicular 
to the plane of the graphene layers and three two-fold 
rotations (C2). There are three irreducible representa- 
tions, two one-dimensional denoted by Tai , and a 
two-dimensional one, Te (see Table IVl)). 



32 (Da) 


E 2C3 3C2 




1 1 1 
11-1 

2 10 



TABLE VI. Character table of the group 32 (D3). 

The basis functions built from the linear combinations 
of Pz and d orbitals which transform as these irreducible 
representations are shown in Table IVIll for ABC trilayer 
and in Table IVllIl for bilayer. 



|*llo), l*2a,-l), 1*2:0,0), l*2:-2,2), l*2:-l,l), l*2;2,-2> 

rA, i*l;;^,o>, 1*5^), 1*1^2.2), l*2;-i.i), 1*2:2,-2) 
r^. {|*?:J.o>,l*?:o,o>}, {l*!I^),l*H!^>} 

{|*^?,),l*?g)}, {|*fl),l*f-l)}, {1*^-2), i*^)} 

{i*^f.),i*i!o)}, {i*ei),i*i!i)}, {1*^,1), i*r-2>} 

{|*M),l*f-l)}> {l*f-l), !*??>}, {1*2^1), !*f-2)} 
{i*f-2),|*fi)} 



TABLE VII. Basis functions for the irreducible representa- 
tions of the small group of the K point for ABC trilayer 
graphene. The basis functions for the K' point can be ob- 
tained by complex-conjugation. 



In the case of bilayer, making then the same steps 
as for trilayer graphene, one finds the symmetry ba- 
sis and intrinsic SOC parameters given in Table IIXI 
Referring to Fig. [Ija), the atomic sites participating in 
the formation of the the split-off bands are denoted by 
Bl and A2, whereas they are labeled as Al and B2 
in Ref. |40, and similarly for the sites contributing to 
the low energy bands. Therefore, to arrive at the SOC 
Hamiltonian shown in Table IIVI (i) one has to rotate 
the symmetry basis into the on-site basis, (ii) re-label 
the sites as Aj ^ Bj, (iii) and make the identifica- 
tion An — — A12 = -Ai/2, Ao = -^{^2/3 — A1/3), 
A4' = — ;^(A2/3 -I- A1/3). The matrix elements are real 
numbers in our calculations because we use different lat- 
tice vectors than in Ref. l40l 



Appendix B 

Here we collect some useful formulae: we give explicitly 
the transformation between the symmetry basis and the 
effective on-site Pz orbital basis for the trilayer case and 
present the low-energy electronic Hamiltonian for bilayer 
graphene. 



1. Transformation between the symmetry basis 
and the on-site effective basis 



The transformation reads: 







/ 


1 


1 











\ 


/ 




\ 








-1 


1 




















l*fff^) 


1 










1 





-1 













"71 










1 





1 


























1 





-1 












V 











1 





1 / 


V 







(Bl) 



Tai 



La, 



i*?:o,o), i*2;l-i>, 1*2:0.0), I*2:i2.2> 



i*to,o), i*2:t-i), i*2:i,o), 1*2:^2,2) 



{|*^o),l*u,)}, {|*^o),l*l!o>}, {l*^,-i), !*!!?>} 

{|*^2),l*f-2)}, {l*f-2),l*i:2)}, {|*2^-l),!*2^})} 
{|*fl),l*f-l>}, {|*f2),|*f-2)} 



TABLE VIII. The same as in Table |VII] but for bilayer 
graphene. 



2. Bilayer graphene low-energy electronic 
Hamiltonian 



For completeness and for comparison with the trilayer 
case given in Eq. ()10|) . we show here the low-energy elec- 
tronic Hamiltonian of bilayer graphene in the on-site ba- 
sis Al, B2. This Hamiltonian has been discuss in many 
publications before, see e.g. the recent review of Ref. |6l|. 
As in Section rVIIi we assume that atoms Al and Bl in 
the layer adjacent to the substrate have a different on-site 
energy A than atoms A2, B2 in the second layer. The 
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I^TaJ - I'^l.O.O/ + 72 «epd-71 I 2.1.-1' ^Epd-Tl I 2.0.0/ 
3 K ■• 



i*r,,,> = l*-i!J) + ^fel*r-i> 



l*r,,,> = l*fo) 



3 VVci,r I g,A2 



^1/2 

^2/3 
^3/3 



9 f . 



3 f 
2 ?d 



5epd(5£pd-7l) 



9 f CVWrT 



TABLE IX. Symmetry basis functions (left column) and intrinsic spin-orbit matrix elements (left column) in terms of Slater- 
Koster parameters for bilayer graphene. We used the notation Vpdcr = V^^^^^ , Sspd — £p — Ed and 71 = 71 + Vdda- 



most important terms arc found to be: 



rVoff 



Hr 



H: 



3uj 



Hy + Ha , 



(vrt)' 



71 V 



7i 







^ 








(: 


:) 







TT 
Tft 



(B2) 
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